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2 I In the framework of /cr-factorization QCD approach and the colour singlet model we 

I consider J/ip inelastic photo- and leptoproduction processes at HERA. We investigate the 

dependences of the single differential and double differential cross section on different forms 
^ '. of the unintegrated gluon distribution. The z and pj. dependences of the spin aligment 

O I parameter a are presented also. Our theoretical predictions agree well with recent data 

^ ' taken by the HI and ZEUS collaborations at HERA. It is shown that experimental study 

of the polarization J/ ijj mesons at low < \ GeV^ is an additional test of BFKL gluon 
Qh! dynamics. 



1 Introduction 



It is known that from heavy quark and quarkonium production processes one can obtain 
unique information on gluon structure function of the proton because of the dominance of 
the photon-gluon or gluon-gluon fusion subprocess in the framework of QCD [1]. Studying 
gluon distributions at modern collider energy (such as HERA, Tevatron) is important for 
prediction of heavy quark and quarkonium production cross sections at future colliders (LHC, 
THERA). At the energies of HERA and LEP/LHC colliders heavy quark and quarkonium 
production processes are so called semihard processes [2-5]. In such processes by definition 
the hard scattering scale ^ ~ mq is large compare to the Aqcd parameter but on the other 
hand /i is much less than the total center-of-mass energy: Aqcd <^ <^ \fs. The last 
condition implies that the processes occur in small x region: x ~ niQ/ y/s <^ 1, and that the 
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cross sections of heavy quark and quarkonium production processes are determined by the 
behavior of gluon distributions in the small x region. 

It is also known that in the small x region the standard parton model (SPM) assump- 
tions about factorization of gluon distribution functions and subprocess cross sections are 
broken because the subprocess cross sections and gluon structure functions depend on a 
gluon transverse momentum kx [2-5] . So calculations of heavy quark production cross sec- 
tions at HERA, Tevatron, LHC and other coUider conditions are necessary to carry out in 
the so called /c^-factorization (or semihard) QCD approach, which is more preferable for the 
small X region than SPM. 

The /cT-factorization QCD approach is based on Balitsky, Fadin, Kuraev, Lipa- 
tov (BFKL) [6] evolution equations. The resummation of the terms In"' [jj? / h^q^Q^y) , 

ln"(^^/AQ(-.j3) ln"(l/a;) and ln"(l/a:) in the /cr-factorization approach leads to the 
unintegrated (dependent from q^) gluon distribution $(a;, q|., //^) which determine the prob- 
ability to find a gluon carrying the longitudinal momentum fraction x and transverse mo- 
mentum qr at probing scale /x^. 

To calculate the cross section of a physical process the unintegrated gluon distributions 
have to be convoluted with off mass shell matrix elements corresponding to the relevant par- 
tonic subprocesses [2-5] . In the off mass shell matrix element the virtual gluon polarization 
tensor is taken in the BFKL form [2-5]: 

L-(g) = (1) 

Nowadays, the significance of the A;T-factorization QCD approach becomes more and 
more commonly recognized [7] . It was already used for the description of a wide class heavy 
quark and quarkonium production processes [8-23]. It is notable that calculations in kx- 
factorization approach provide results which are absent in other approaches, such as the fast 
growth of total cross sections in comparison with SPM, a broadening of the pt spectra due 
to extra the transverse momentum of the colliding partons and other polarization properties 
of final particles in comparison with SPM. 

We point out that heavy quark and quarkonium cross section calculations within the SPM 
in the fixed order of pQCD have some problems. For example, the very large discrepancy (by 
more than an order of magnitude) [24, 25] between the pQCD predictions for hadroproduc- 
tion J / ip and T mesons and experimental data at Tevatron was found. This fact has resulted 
in intensive theoretical investigations of such processes. In particular, it was required to use 
additional transition mechanism from cc-pair to the J/ if) mesons, so-called the colour octet 
(CO) model [26], where cc-pair is produced in the color octet state and transforms into 
final colour singlet (CS) state by help soft gluon radiation. The CO model was supposed 
to be applicable to heavy quarkonium hadro- and leptoproduction processes. However, the 
contributions from the CO mechanism to the J/ip meson photoproduction contradict the 
experimental data at HERA for z distribution [27-30] . 

Another difficulty of the CO model arc the J /ip polarization properties in pp-interactions 
at the Tevatron. In the framework of the CO model, the J/ijj mesons should be transverse 
polarized at the large transverse momenta p^. However, this is in contradiction with the 
experimental data, too. 

The CO model has been applied earher [31, 32] in an analysis of the J/ip inelastic produc- 
tion experimental data at HERA [33] . However, the results do not agree with each other [32] . 
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We note that the shapes of the Q^, and y* distributions are not reproduced by the calcu- 
lation [31], and the z distributions [32] contradict the HERA experimental data too. Results 
obtained within the usual coUinear approach and CS model [34-37] underestimate experi- 
mental data by factor about 2. 

The inelastic J/ip production at HERA in the CS model with fcT^-factorization also was 
considered in [15, 16]. The results [16] agree with HI experimental data [33] both in normal- 
ization and shape only at quite small charmed quark mass rric — 1.4 GeV. The theoretical 
prediction [15] are stimulated the experimental analysis of J/%1) polarization properties at 
HERA conditions. 

Recently the new experimental data on the inelastic J/^ photo- and leptoproduction 
at HERA were obtained by the HI [38, 39] and ZEUS [40] collaborations with increased 
statistics and precision as compared with previous experimental analyses [33] . Based on the 
above mentioned results here we will use the CS model and the fc^^-factorization approach for 
the analysis of the data [38-40]. We investigate the dependences of the single differential and 
double differential J / ip production cross section on different forms of the unintegrated gluon 
distribution. Special attention is drawn to the unintegrated gluon distributions obtained 
from BFKL evolution equation which has been applied earlier in our previous papers [12- 
16]. For studying J/ip meson polarization properties we calculate the z and pr dependences 
of the spin aligment parameter a. 

The outline of this paper is as follows. In Section 2 we present, in analytic form, the total 
and differential cross section for the inelastic J/ip photo- and leptoproduction in the CS 
model with /c^-factorization, and give the formulas for the relevant partonic subprocess off 
mass shell matrix elements. In Section 3 we present the numerical results of our calculations 
and compare them with the HI [38, 39] and ZEUS [40] data. Finally, in Section 4, we give 
some conclusions. 



2 Analytic results 

In this section we calculate total and differential cross section for inelastic J/ijj photo- 
and leptoproduction in the CS model with /cr-factorization, and give the formulas for the 
relevant partonic subprocess off mass shell matrix elements. 

2.1 Kinematics 

As indicated in Fig. 1, we denote the 4-momenta of the incoming electron and proton 
and the outgoing electron, proton remnant, J/ip meson and gluon by p^, Pp, p^, p'p, Pip and 
Pg, respectively. The initial virtual photon and BFKL gluon have a 4-momenta qi = Pe—p'e 
and q2 = Pp — p'p, so that 4-momentum transfer Q"^ = —qf. In our analysis below we will use 
the Sudakov decomposition, which has the following form: 

Py, = aiPe + PlPp + P^T: Pg = 0£2Pe + (^2Pp + PgT, 

qi = XiPe + qiT, ?2 = X2Pp + q2T, 

where p^r, Pgr, Qit and q2T are transverse 4-momenta of corresponding particles, and 

2 2 2ri22 22 /o\ 

= "^V" Pg = 0' ^1 = llT: ?2 = ?2T- (3) 
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In the ep cm. frame we can write: 

Pe = ^rs|'2 (1, 0, 0, 1), pp^ (1, 0, 0, -1), (4) 

where we neglect the masses of the electron and proton. The Sudakov variables are expressed 
as follows: 

Oil = exp(y^), a2 = ^-^^ exp(yg), 

m II 

where m'^j, = + P^T) Vtp are the rapidities of J/%1) meson and final gluon respec- 

tively in the ep cm. frame. Prom conservation laws we can easy obtain following conditions: 

Xl=ai + a2, X2=l3i+I32, (\iT + ^.2T = Vi>T + VgT- (6) 

Also the variable z = {p^p ■ pp) / {qi ■ Pp) is used for a description of quarkonium photo- and 
leptoproduction processes. In the rest frame of the proton one has z = E^/Ej. 

2.2 Inelastic J/ip leptoproduction cross section 

In the /cy-factorization approach the differential cross section for inelastic J/ip leptopro- 
duction may be written as: 

da{ep ^ e' J/iP X) ^ ^ <^{x2, qlr, ^ dqjj, da{e g* ^ e' J/V' g'), (7) 

X2 ^TT 

where 02 is initial BFKL gluon azimuthal angle, $(a;2, qij^, A*^) is an unintegrated gluon 
distribution in the proton. The eg* ^ e' J /if) g' cross section is given by: 

da{eg* ^ e'J/i^g') = ^ Y.W\lnK{eg* - e'J/^l^g') x 

^X2S 

C^Ve d^PiJ d?Pg r(4)/ , / N 

^ (2^0^ j2nfW, ' + - - - ' 

where I^Isha(^5'* ~^ ^' J/'4^g') is the off mass shell matrix element. In (8) indicates 
an averaging over and a sum over the final polarization states. Prom (7) and (8) we obtain 
the following formula for the inelastic J/-^ leptoproduction differential cross section in the 
/cr-factorization approach: 

/ ^; , 1 ^{x2, c&T, 1^'^) dz 

da{ep ^ e' J ijj X) = -— , ^ ' ^ \ — dy^ x 

USn-^ {X2 sy {1 - xi) z (1 - z) 

X E mluAieg* - e' J/^ g') dpj^ dQ^ dqj^ ^ ^' 

where (pi and 0^ are azimuthal angles of the initial virtual photon and J/ip meson respec- 
tively. 
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2.3 Inelastic J/'ijj photoproduction cross section 

As in leptoproduction case, in tlie /cr-factorization approach the differential cross section 
for inelastic J/il) photoproduction may be written as: 

d<7(7P- J/V^) = — $(X2, q^T, ^i')^dql^da{^g* ^ J/i^g'). (10) 

If we take the limit — > and Xi — 1, we easy obtain the following formula for the inelastic 
J/ip photoproduction differential cross section in the /cT-factorization approach by analogy 
with the leptoproduction case: 

1 dz 

We note that formulas for the the differential cross section for inelastic J/ijj photo- and 
leptoproduction in the usual parton model may be obtained from (9) and (11), if we take 
the limit q^j^ — > and average them over the transverse momentum vector q2T- 

2.4 Off mass shell matrix element 

There are six Feynman diagrams (Fig. 2) which describe partonic subprocess 7(7* —> 
J/ipg' at leading order in as and a. In the framework of the CS model and the nonrela- 
tivistic approximation the production of the J/ip meson is considered as a production of a 
quark-antiquark system in the colour singlet state with orbital momentum L — and spin 

momentum S = 1. The binding energy and relative momentum of the quarks in the J/i/j 
meson are neglected, resulting in = 2mc, where is charm mass. The amplitude of 
the process j g* — > J/ip g' may be obtained from the amplitude of the process 75'* ccg' 
after replacement: 

v{pc) u{pc) J {p^) = e(p^) (p^ + m^) (12) 

2Jm^ w3 



where Pc = Pxi,/^, e{p^) is a 4-vector of the J /ip polarization, 1/a/3 is the color factor, ■0(0) 
is the nonrelativistic meson wave function at the origin. The matrix element is: 



xSp 



M ^Ccg^ e^(gfi) ea{q2) ep{pg) x 



(Pc - gi)^ - ml {-pc - PgY - m 



+ 5 permutations of all gauge bosons. Here ^^{qi) and e^{q2) are polarization vectors of 
the initial photon and gluon respectively, e^ipg) is a 4-vector of the final gluon polarization. 
The summation on the meson and final gluon polarizations is carried out by covariant 
formulas: 
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The initial BFKL gluon polarization tensor is taken in form (1). For the photon we use the 
usual expression 

T.^'i^^y''i^i)--9^'' (16) 

in photoproduction case and the full lepton tensor (including also the photon propagator 
factor and photon-lepton coupling) in leptoproduction case: 

T.e'{qi)e*^{qi) = 2^ (^-g^^ + ^^ . (17) 

For studying J/ip polarized production we introduce the 4- vector of the longitudinal 
polarization ei{p^) as follows [41]: 

The evaluation of X! |Af Hha for photo- and leptoproduction cases was done analytically 
by the REDUCE program. Also in our calculations wc have used the JB [42] and KMS [43] 
parametrizations of the unintegrated gluon distributions (see also [7] for the detail informa- 
tion) . 



3 Numerical results 

In this section we present the theoretical results in comparison with recent experimental 
data taken by the HI [38, 39] and ZEUS [40] coUaborations at HERA. 

There are three parameters which determine the common normalization factor of the cross 
section under consideration: J/i(j meson wave function at the origin ijj{0), charmed quark 
mass rric and factorization scale fi. The value of the J/ip meson wave function at the origin 
may be calculated in a potential model or obtained from the well known experimental decay 
width r(J/^ fi+ In our calculation we used |^(0)p = 0.0876 GeV^ as in [44]. 

Concerning a charmed quark mass, the situation is not clear: on the one hand, in the 
nonrelativistic approximation one has = m,^/2 = 1.55 CeV, but on the other hand there 
arc examples when smaller value of a charm mass nic = 1.4 GeV is used [32, 45]. However, 
in our previous paper [16] we analyzed in detail the influence of charm quark mass on 
the theoretical results. We found that the main effect of change of the charm quark mass 
connects with final phase space of J/ijj meson, and in the subprocess matrix elements this 
effect is neglectable. Taking into account that the value of rric = 1.4 GeV corresponds to the 
unphysical phase space of J/ip state, in the present paper we will use value of a charm mass 
rric — 1.55 GeV only. 

Also the most significant theoretical uncertanties come from the choice of the factorization 
scale fip and renormalization one yU^. One of them is related to the evolution of the gluon 
distributions ^{x, Qt; A^f)? the other is responsible for strong coupling constant as{iJ>%). As 
often done in literature, we set /ip = /ir = fi. In the present paper we used the following 
choice //^ = q|j, as in [16, 46]. 
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3.1 Inelastic J/^ leptoproduction at HERA 

The integration limits in (9) are taken as given by kinematical conditions of the HI 
experimental data [39]. One kinematical region^ is 2 < < 100 GeV^, 50 < ly < 225 GeV, 
0.3 < z < 0.9, p^?p > IGeV^ and other kinematical region is 12 < < 100 GeV^, 
50 <W < 180 GeV, p^^ > 6.4 GeV^ 0.3 < 2 < 0.9 and p;?p > 1 GeV^ Here and in the 
following, we used Aqcd = 250 MeV. 

The results of our calculations are shown in Fig. 3 — 5. Fig. 3 shows the single differential 
cross sections of the inelastic J/ip meson leptoproduction obtained in the first kinematical 
region at ^ = 314 GeV. Curve 1 corresponds to the SPM calculations at the leading order 
approximation with the GRV (LO) gluon density, curves 2 and 3 correspond to the kx- 
factorization results with the JB (at A = 0.35 [17, 23]) and the KMS unintegrated gluon 
distributions. One can see that results obtained in the CS model with /cT-factorization agree 
very well with the HI experimental data. The SPM calculation are lower than the data by 
a factor 2 — 3. 

We would like to note the difference in the transverse momenta distribution shapes be- 
tween curves obtained using the fcr-factorization approach and the SPM. This difference 
manifests the pt broadening effect which mentioned earlier. It is visible also that only the 
fcr-factorization approach gives a correct description of the p^j- spectra. However, we note 
that the p^^ distributions somewhat less well described (in contrast with p^j, spectra) at 
the large values of the J/ip transverse momenta (see Fig. 3d). 

Also we point out the good description of the z distributions which obtained in the /ct- 
factorization approach in contrast with CO model results [32], except for the region z < 0.3, 
where the contribution of the resolved photon process may be large [47]. 

Fig. 4 shows the single differential cross sections of the inelastic J/ip meson production 
obtained in the second kinematical region at y/s = 314 GeV. Curves 1 — 3 are the same as 
in Fig. 3. We find also good agreement between results obtained in the CS model with kr- 
factorization and HI data. It is notable that in this kinematical region in contrast with first 
one the both p^-^ and p^?^ transverse momenta distributions agree well with the experimental 
data. 

The double differential cross sections da/dQ'^dz and da/d'p*^rpdz (Fig. 5) obtained with 
/cT-factorization in the different z regions 0.3 < z < 0.6 (Fig. 5a, b), 0.6 < 2; < 0.75 
(Fig. 5c, d) and 0.75 < z < 0.9 (Fig. 5e, f) agree with the HI data. We note that double 
differential cross sections rfcr/rfp^^ dz somewhat less well described at the large z (see Fig. 5f). 
However, in this region the contribution of the diffractive processes may be large. All of these 
contributions are not in our consideration. 

It is interesting to note that results obtained with the JB unintegrated gluon distribution 
at A = 0.35 and the KMS ones, which effectively included about 70% of the full NLO 
corrections to the value of A [43], coincide practically in a wide kinematical region. 

Fig. 3 — 5 show that the /cT-factorization results for inelastic J /ip leptoproduction with 
realistic value of a charm mass ttlc = 1.55 GeV agree well with the HI experimental data 
without any additional cc — J/ip fragmentation mechanisms, such as the CO contributions. 

3.2 Inelastic J/ip photoproduction at HERA 

The integration limits in (11) are taken as given by kinematical conditions of the HI [38] 

•^Here we denote the J/i/i meson transverse momentum and rapidity in the 7* p cm. frame by p'^j, and j/J, respectively. 
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and the ZEUS [40] data. One kinematical region which corresponds to the HI experiment is 
60 < < 240 GeV, 0.3 < z < 0.9, 1 < p^,. < 60GeV^ and other kinematical region which 
corresponds to the ZEUS experiment is 50 < < 180 GeV, OA < z < 0.9, p^^ > 1 GeV^ 

The resuhs of our calculations are shown in Fig. 6 — 8. Fig. 6 and 7 show the total and 
single differential cross sections of inelastic J/ip meson photoproduction in comparison with 
the HI and ZEUS data, respectively. As in previous section, curve 1 corresponds to the SPM 
calculations at the leading order approximation with the GRV (LO) gluon density, curves 
2 and 3 correspond to the /ct- factorization results with the JB (at A = 0.35 [17, 23]) and 
KMS unintegrated gluon distributions. 

The W dependences of the total J/ip photoproduction cross section at 0.3 < z < 0.9, 
0.3 < z < 0.8 and 0.4 < z < 0.9 are plotted in Fig. 6a, Fig. 6b and Fig. 7a respectively. One 
can see that results obtained in the CS model with /c^-factorization agree very well with the 
HI [38] and the ZEUS [40] experimental data. The SPM results are lower than the data by 
a factor 2. 

Concerning the shapes of the p^y distribution (Fig. 6c and 7b), one can note a difference 
between the /cT-factorization and the SPM curves. As in leptoproduction case, this difference 
manifests the pr broadening effect which was mentioned earlier. It is visible also that only 
the /cr-factorization approach gives a correct description of the HI data. 

The z distributions are shown in Fig. 6d, e, f and Fig. 7c, d, e at different p^r cuts in 
comparison with the HI and ZEUS data, respectively. One can see that good agreement be- 
tween the /cT-factorization curves and the experimental data is observed. The z distribution 
somewhat are less well described at p^T > 3 GeV (see Fig. 6f). The dispcrancc between the 
leading order SPM calculations and the experimental data is about factor 2 at p^j- > 1 GeV 
and about order of magnitude at p^T > 3 GeV for z ~ 0.8. Also we note that in the region 
2; < 0.3 the contribution of resolved photon process may be large [47], as in leptoproduction 
case. 

The double differential cross sections da/dp'^j^dz (Fig. 8) in the different z regions 0.3 < 
z < 0.6 (Fig. 8a), 0.6 < ^ < 0.75 (Fig. 8b) and 0.75 < z < 0.9 (Fig. 8c) are well described 
by the /cr-factorization approach. 

It can be seen that the results obtained with the JB unintegrated gluon distribution 
with A = 0.35 and the KMS ones (which effectively included the main part of the full 
NLO corrections to the value of A) practically coincide in a wide kinematical region, as in 
leptoproduction case. 

Fig. 6 — 8 show that the /ct- factorization results for inelastic J/t/j photoproduction with 
realistic value of a charm mass rric = 1.55 GeV agree well with the HI and ZEUS experimental 
data without any additional cc ^ J/ip fragmentation mechanisms, such as CO contributions. 



3.3 Polarization properties of the J/ip meson at HERA 

As it mentioned above, one of differences between the /c^-factorization approach and the 
SPM is connected with polarization properties of the final particles. In the present paper 
for studying J/ip meson polarization properties we calculate the p^ and z dependences of 
the spin aligment parameter a [14-16]: 

^ ^ da/dw + daL/du ^ ' 
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where ul is the production cross section for the longitudinally polarized J/il) mesons, u) 
p^T, z. The parameter a controls the angular distribution for leptons in the decay J/ip - 
/J,'^ /I" (in the J/ip meson rest frame): 

dr(j/^ ^ 2 



dcos9 



~ 1 + a cos^ 9. (20) 



The cases a — 1 and a — —1 correspond to transverse and longitudinal polarization of the 
J/ip meson, respectively. 

In our previous paper [16] we analyzed in detail the and p^^. dependences of the 
spin parameter a in leptoproduction case. We found that it is impossible to make of exact 
conclusions about a BFKL gluon contribution to the polarized J/ ip production cross section 
because of large additional contribution from initial longitudinal polarization of virtual pho- 
tons. However at low and in photoproduction limit these contributions are negligible. 
This fact should result in observable spin effects of final J/ip mesons, connected with the 
/cT-factorization effects. In this paper we have performed such calculations for the inelastic 
J/ip photoproduction process. 

The results of our calculations are shown in Fig. 9 and 10. Fig. 9 shows the parameter 
q; as a function z and p^T in comparison with the HI experimental data which obtained in 
the kinematical region 60 < W < 240 GeV, 0.3 < ^ < 0.9 and 1 < p^^ < 60 GeVl Curve 
1 corresponds to the SPM calculations in the leading order approximation with the GRV 
(LO) gluon density, curve 2 corresponds to the /cT^-factorization results obtained with the JB 
(at A = 0.35 [17, 23]) unintegrated gluon distribution. One can see that the z dependence 
of the spin parameter a is not sensitive to the results of different approaches, included the 
nonrelativistic QCD predictions (see [38]). However, the behavior of the q;(pt) is different in 
the A;T-factorization approach and the SPM (see Fig. 9b). Although the experimental points 
have large errors they tends to support the fc^-factorization theoretical predictions. 

Fig. 10 shows the p^T dependence of the spin parameter a in comparison with the ZEUS 
experimental data which obtained in the kinematical region 50 < W < 180 GeV, 0.4 < z < 
0.9 (Fig. 10a and 10c) and OA < z < 1 (Fig. 10b and lOd). We note that in Fig. 10a and 
10b the quantisation axis is chosen to be opposite of the incoming proton direction in the 
J/ip rest frame, 6 is the opening angle between the quantisation axis and the fi^ direction of 
flight in the J/ijj rest frame. This frame is known as the "target frame" [40]. In Fig. 10c and 
lOd, the quantisation axis was defined as the J/t/j direction of fiight in the ZEUS coordinate 
system. This frame is known as the "hehcity basis" [40, 48]. Curves 1 and 2 are the same 
as in Fig. 9. 

It is visible that only the fcr-factorization approach gives a correct description of the 
ZEUS data, although the experimental points have large errors. We also have large difference 
between predictions of the leading order of SPM and the /cr-factorization approach. The 
SPM predictions lies somewhat below the data at low p^T and somewhat above at high 
p^T- Therefore experimental measurement of polarization properties of the J/ip mesons will 
be an additional test of BFKL gluon dynamics. 



4 Conclusions 

In this paper we considered the inelastic J/ip meson photo- and leptoproduction at HERA 
in the colour singlet model using the standard parton model in leading order in as and the 
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/cT-factorization QCD approach. We investigated the total cross section, single differential 
and double differential cross sections of inelastic J /if) production on different forms of the 
unintegrated gluon distribution. The and z dependences of the spin aligment parameter a 
presented also. We compared the theoretical results with recent experimental data taken by 
the HI and ZEUS collaboration at HERA. We have found that the /c^-factorization results 
(in contrast with the SPM ones) with the JB and KMS unintegrated gluon distributions 
agree well with the experimental data at realistic value of a charm mass mc = 1.55 GeV, 
1^(0) P — 0.0876 GeV^ and Aqcd = 250 MeV without any additional transition mechanism 
from cc-pair to the J/ip mesons (such as given by the CO model). We also found that 
results obtained with the JB unintegrated gluon density at A = 0.35 and KMS one, which 
effectively included about 70% of the full NLO corrections to the Pomeron intercept A, 
practically coincide in a wide kinematical region for J/ip production processes at HERA 
conditions. Finally, it is shown that experimental study of a polarization of J/ip meson at 
low < 1 GeV^ should be additional test of BFKL gluon dynamics. 
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Fig. 1. Diagram for ep —>■ e' J/ipX process. 

Fig. 2. Feynmans diagram used for description partonic process 75' —> J/i^ g' process. 

Fig. 3. The single differential cross sections of the inelastic J/il) leptoproduction obtained 
in the kincmatical region 2 < < 100 GeV^ 50 < < 225 GeV, 0.3 < z < 0.9 and 
p^^ > 1 GeV^ at ^/s = 314 GeV in comparison with the HI [39] data. Curve 1 corresponds 
to the SPM calculations at the leading order approximation with GRV (LO) gluon density, 
curves 2 and 3 correspond to the /cT-factorization QCD calculations with JB and KMS 
unintegrated gluon distribution. 

Fig. 4. The single differential cross sections of the inelastic J/ijj leptoproduction obtained 
in the kinematical region 12 < < 100GeV^ 50 < < 225 GeV, p^^ > 6.4GeV^ 
0.3 < ^ < 0.9 and p*^\. > 1 GeV^ at = 314 GeV in comparison with the HI [39] data. 
Curves 1 — 3 are the same as in Fig. 3. 

Fig. 5. The double differential cross sections of the inelastic J/^jJ leptoproduction obtained 
in the kinematical region 2 < < 100 GcV^ 50 < W < 225 GeV, 0.3 < ^ < 0.9 and 
p^|. > 1 GeV^ at ^/s — 314 GeV in comparison with the HI [39] data. Curves 1 — 3 are the 
same as in Fig. 3. 

Fig. 6. The total and single differential cross sections of the inelastic J/ip photoproduction 
obtained in the kincmatical region 60 <W < 240 GeV, 1 < p^^ < 60 GeV^ 0.3 < ^ < 0.9 
in comparison with the HI [38] data. Curves 1 — 3 are the same as in Fig. 3. 

Fig. 7. The total and single differential cross sections of the inelastic J/ ip photoproduction 
obtained in the kinematical region 50 < < 180 GeV, p^^, > 1 GeV^, 0.4 < ;2 < 0.9 in 
comparison with the ZEUS [40] data. Curves 1 — 3 are the same as in Fig. 3. 

Fig. 8. The double differential cross sections of the inelastic J/ip photoproduction obtained 
in the kinematical region 60 < W < 240 GeV, 1 < p^^ < 60 GeV^ 0.3 < ^ < 0.9 in 
comparison with the HI [38] data. Curves 1 — 3 are the same as in Fig. 3. 

Fig. 9. The parameter a as a function z and p^^ for the inelastic J/ip photoproduction 
process which obtained in the kinematical region 60 < W < 240 GeV, 0.3 < ^ < 0.9 
and 1 < p^y < 60GeV^ in comparison with the HI [38] data. Curve 1 corresponds to 
the SPM calculations at the leading order approximation with GRV (LO) ghion density, 
curve 2 corresponds to the /cT-factorization QCD calculations with JB unintegrated gluon 



distribution. 



Fig. 10. The parameter a as a function p^^ for the inelastic J photoproduction process 
which obtained in the kinematical region 50 < < 180 GeV, 0.4 < z < 0.9 (Fig. 10a, c), 
0.4 < 2; < 1 (Fig. 10b, d) and 1 < p^^, < 60 GeV^ in comparison with the ZEUS [40] data. 
Curves 1 and 2 are the same as in Fig. 9. 
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